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The purpose of the study undertaken as part of this grant is to conduct 
research for a better understanding of the mechanism of heat transfer in the 
lunar regolith. We worked on a method to make measurements of the thermal 
conductivity of lunar core samples in the laboratory. In the period covered 
by this proposal, we concentrated our effort on the experimental design and 
the development of a theory that is necessary to interpret the experimental 
data. 

After examining several possible methods of measurement, we 
concluded that one of the best methods that can be applied to the measurement 
of thermal conductivity of lunar core samples is to make the measurements 
while the lunar material is still in the core tube. This approach would reduce 
the possibility of physical and chemical disturbances to the sample. The 
approach we chose is to heat the sample externally by radiation at a known 
rate, measure the variation of temperature at the surface of the core sample, 
and determine thermal conductivity by comparing the observed temperature 
with the theoretically expected one. 

We chose to heat the sample by radiation because a direct contact of the 
heating element to the sample can be avoided. The thermal conductivity of 
lunar regolith material is very sensitive to the interstitial gaseous pressure. 

To simulate the lunar surface condition, the measurement must be made at a 

-3 -3 

pressure below 10 torr. At gaseous pressures below 10 torr, the heat 

transfer due to convective motion of the gas is negligible, so that the rate of 
thermal energy transmitted from the heating unit to the samplt can be esti- 
mated accurately by the theory of thermal radiation. 

One of the most important tasks of this project was to obtain a mathe- 
matical expression for a surface temperature of core sample changes with 
time when it is heated externally at a known rate. The sample of lunar 
regolith material is contained in a metallic circular cylinder. Therefore, 
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the core sample constitutes concentric composite cylinders of two 
materials. The variation of temperature in this system, subject to the 
inflow of heat through the outer surface of the cylinders, should be given. 
The solution will be a function of the physical properties of the materials 
consisting of the outer and the inner cylinders. Since the physical pro- 
perties of the metallic core tube are known, those of the lunar regolith 
material will be obtained by comparing the solution with experimentally 
observed temperature. 

As far as we know, the mathematical solution to this problem 
cannot be found in the literature. We spent a considerable amount of 
time to formalize the problem, obtain a solution, and put it in a form 
readily applicable to the analysis of data. The details are given in 
Appendix J. Computer programs were constructed so that the least- 
square fit to the data is done by trial and error to yield a best estimate 
of the thermal conductivity of the lunar regolith sample. 

Two types of lunar core samples are available for our experiment. 
One is the Apollo 15 type drive-tube core sample. The other is the Apollo 
15 type drill-core sample. To provide radiation to the outer wall of the 
core sample, two sets of a heating unit were constructed to fit each type 
of core sample. The heating units are coaxial cylinders. Heat transfer 
from the heating unit to the core sample is two-dimensional thermal 
radiation. The heat flux Q, per unit time, per unit axial length of the 
core tube is 


Q = A 

c 




) 


( 1 ) 
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where T and T are the surface temperatures o£ the core tube and the 
c h 

heater, * = (<r /a.) and €, = (<r. /tr. ) are the thermal emissivities of 

c c b h h b 

the core tube and the heater surfaces, A is the area of the core tube 

c 

surface per unit length of the core tube, J is the ratio of the surface areas 

of the core tube and the heater per unit axial length of the cylindrical system. 

If the radius of the core tube is r and the inner radius of the heater r, , 

c h 

then 


A = 2 Tr 
c c 


and 


6 - Z tr r /2 It r. = r / r, 
c h c h 

When the temperature difference between the heater and the core 
tube is small, (1) is approximated by 


4A T 

c c 

— — 

a r a, 
c h 


-AT 


( 2 ) 


where 


AT 



- T 


c. 


Since «r, and a are constants that depend on the surface finish of the 
h c 

materials in a complicated way, the coefficient in (2) 


f = 


4A T 


3 


c c 




(3) 


can be determined experimentally. A co r e tube, not filled with sample 

material, is used for this experiment. Let the specific heat, the density 

the outer and the inner radii of the core tube bee,, P ,, r ( = r ) and r . 

2 2 2 C I 

Then, the heat capacity of the core tube per unit length of the core tube is 
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2 2)c 0 _ . Since the core tube is made of thermally conductive 

▼ (r_ - r 2 ~2 

metai and the rate of temperature variation is comparatively small, it can 
be assumed that the temperature of the core tube is uniform and equal to 
its surface temperature. Therefore, the rate of heat inflow into the core 
tube is given by 


dT 

Q = IT (r - r ^)c s * 
'2 1 ’2 2 dt 


(4) 


which is to be equated to the rate of heat transfer from the heater to the 
core tube. Equation (4) combined with (2) allows us to calculate f. The 
experimental result showed that a better fit to the data is obtained if an 
additional term, proportional to the core tube temperature, is included in 
(2), i. e. , 


Q = f AT + g T 

c 


(51 


Equating (4) and (5) and integrating with respect to time, we obtain 

7T (r, 2 - r, 2 )c « T (t) = f ( AT (t) dt + g( T (t) dt (6) 


Values of f and g are determined from the measured T and T by the 

c n 

least square criteria. 

Comparison of the experimentally determined f with the theoretical 

expression (3) yields an estimate of « . The surface of the heater is 

c 

painted a non- glossy black. Therefore, 4 =1. The estimated t is 

h c 

0.47 for the Apollo 15 type drive-tube and is 0. 64 for the Apollo 15 type drill 
core tube. The material of the Apollo 15 type drive-tube is an aluminum 
alloy 6016-T6 with a surface finish of sulfuric anodize type II, class 1 per 
mil - A - 8625 manufactured by J. T. Ryerson & Co. The surface layer is 
less than 0. 001" thick and consists mainly of aluminum oxide. The 



The Apollo 15 type drill core-tube is made of a titanium alloy Ti-6A1-4V. 
The estimated emissivities are in reasonable agreement with published 
emissivitles for these materials. 

The experimentally determined g has a negative sign. This 
implies that the second term on the right-hand side of (5) represents a 
loss of heat from the core tube, probably by conduction along the axial 
direction of the core tut?. 

Experimental Setup 

Figure 1 illustrates the heater-sample holder assembly for the 
experiment. The heater for the measurement of the Apollo 15-type drive- 
tube core sample is 8. O' 1 long, 2. 135" in inner diameter and 2. 677" in 
outer diameter. The heater for the measurement of the Apollo 15-typ** 
drill-core sample is 8. 0" long, 1. 40" in inner diameter and 1. 937" in 
outer diameter. The sample holders are fabricated in accordance with the 
sizes of the core tube, 14. 75" long and 4. 39 cm in outer diameter for the 
Apollo 15-type drive-tube and 16.75" long and 2. 33 cm in outer diameter 
for the Apollo 15-type drill-core. Since the heater and the sample holder 
will be in contact with the lunar samples, materials permissible by NASA 
standards are used for the fabrication of the assembly. For temperature 
measurements on the surface of the core tube, thermocouples sheathed in a 
copper disk, 5 mm in diameter and having the same curvature as the outer 
surface of the core tube, were *ied on the core tube by a polyimide strap 
and springs. For the temperature measurement at the surface of the 
heater, stainless steel thermocouples, sheathed in a copper strip, were 
bolted to the inner surface of the heater. 

The heater- sample holder assembly is encased in a vacuum 
chamber 36. 0" long and 5. 90" in inner diameter. There are two reasons 
for making the experiment in a vacuum. First, since the thermal conduc- 
tivity of the lunar regolith material strongly depends on the interstitial 
gaseous pressure, gas pressure in the sample must be kept below 10**(-3) 
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torr during the measurement to simulate the lunar, surface condition. 

Second, to simplify the theoretical interpretation of heat transfer from 
the heater to the core tube, heat transfer due to convective motion of 
the gas filling the space between the heater and the core tube must be 
eliminated by keeping the gas pressure as low as possible. 

Our vacuum system is schematically shown i.. Figure 2. It is 
necessary to control the time rate of evacuation so that a reduction in 
pressure does not cause any mechanical disturbance in the sample. A 
fine metering valve, illustrated in Figure 2, is used to establish and 
control the very slow rate of evacuation. An optimum pumping rate was 
sought experimentally. A glass tube of the same dimension as the core 
tube was filled with material simulating lunar soil and artificial stratifications 
were constructed by intercalating layers of clay of a different color than 
the simulated lunar material. The adoption of a rate of 760 torr/12 hours 
resulted in no disruption in the strata and no leakage of soil particles 
through the small hole in the plug inserted in the top end of the core tube. 

Test Measurement 

A series of test experiments was performed on simulated samples 
of lunar regolith material. The result of an experiment with glass beads 
and powdered Knippa basalt, in an Apollo 15-type drive tube, is summarized 
in Horai et al. (1974)!*' An experiment with Apollo 12 simulant lunar regolith 
material, contained in an Apollo 15-type drill-core tube, also yielded a 
satisfactory result. We think that the feasibility of the method developed 
under this proposal has been demonstrated. 

As a continuation of this study, we plan to make measurements on 
lunar core samples. We would like to note that the Lunar Sample Analysis 
Planning Team has already approved the usage of one Apollo 15-type drive- 
tube core sample and one or two Apollo 15-type drill-core samples 
returned by Apollo 16 and 17 missions. The performance of the measure- 



ments and the interpretation of the results will be undertaken under the 
companion contract (NGR33-008- 1 69). 

The work performed under this grant is being continued under a 
renewal of grant NGR33-008- 177. 


* Horai, M. Langseth, Jr. , A. Wechsler, J. Winkler, D. Colvin and S. 

Keihm; (1974) A new technique of thermal conductivity measure- 
ment of lunar core samples; The Moon, v. 9, p 243. 
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Ue want to dlacuaa a conduction of heat in coopoaite coaxial circular cylinders. 


In cylindrical coordinate system, the inner cylinder (0 * r 6 r^) is of on» substance 
with thermal conductivity k^, specific heat c^, density and thermal diffusivitj 
(= - the outer cylinder (r, - r - r ? ) is of another cubetmce with thermal con- 

n. 1 1 k- 

ductivity k 2 , specific heat c 2 , density f 2 and thermal diffuaivity * 2 ( s — - — — ). The 
cylinders extend to infinity in both positive and negative z-dir«ctions. Initially, the 
cylinders are at constant temperature zero and for time t > 0 heated externally at a 
constant rate Fq per unit time and per unit length of the cylinders. 

In the case we are going to discuss, the outer cylinder (metal) is 10^ times more 
conductive than the inner cylinder (lunar regolith material under high vacuum). Besides, 
we are going to adopt a comparatively small "ate of heating Fp. We will assume that the 
outer cylinder is a perfect thermal conductor (k 2 *<*>) and no temperature difference 
exists between the outer and the inner surfaces of the outer cylinder. 

Therefore, the equation to be solved for the temperature v^(r,t) that is defined 
in the ranges of parameters 0 - r - r^ and t > 0 is 


1ZL 

?r‘ 


/ 

r 


Jr 


J v, 




- o 


( 1 ) 


with the initial and the boundary conditions 


v, (r,0) = 0 


( 2 ) 


and 


K * r 2 “ r l * c 2 r 2 




+ 2£r i k i — 


■r, 


ft, 


= F, 


( 3 ) 


The temperature in the outer cylinder v 2 (r,t), defined in the ranges of parameters 
- r - r- and t > 0, is constant as a function of r and is equal to v^(r^,t). 


A solution of (1) with the conditions (2) and (3) will be obtained by the 
Laplace transformation method (see, for example, Carslaw and Jaeger (1959, chapter 12)). 
If we denote the Laplace transformation of *^(r,t) by 

^(r.p) =^exp(-pt) v^(r,t) dt 

,the subsidiary equation that is obtained by the Laplace transformation of (1) becomes 

i dv. 


♦ 


dr 


r dr 


** - \\ * 0 


(4) 


where 


9 X = p/* x 

with the boundary condition 

* 

K(u'-rS)c*f*f*. 


♦ -2 a. y , 


A* 






- F»/f 


(5) 


The solution of (4) that satisfies the condition (5) is 


v x (r,p) 





&(*•- r,‘) f * c, y> j t » Jzr.A, 1, (f,n) 


(6) 


where I Q (z) and I^(z) are the zero order and the first order modified Bessel functions. 

The solution of (1) with the conditions (2) and (3) is obtained by the inversion 
theorem for the Laplace transformation 




v 1 (r,t) = 


2~i 




( r, > ) A a 


r 


(7) 


where v^r,/*) is the function obtained from v^r,?) by replacing real variable p by 
complex variable A . 




Since the integrand of (7) haa no branch point at A =0, the integral can be 

evaluated by calculating the reaiduea at the polea of the integrand. Inapection of (6) 

2 

reveala that the polea exlat at A = 0 and A = - where ia the n'th poaitive 

root of 

, (r 2 2 - r. 2 ) C 

i j* ' * (*! r„) J o< r l + J l< r l * 0 («> 

^ ^1 6 1 

Here, Jq(x) and J^(x) are the zero order and the first order Beaael functions. 

If we put = £r^ f i c 2 and D 2 = ^ ( r 2^” r l^ ? 2 C 2 * the resldual at A = 0 




t / t'p.c. 

* Pi ** At p t t Pi 


/ _ p, ( p> » ■» Pi) ) 

'**' (p.’P.Y ) 


and that at A s - ^ / 

± n 




From (9) and (10) we have 


7 - f £ j Kt % M> _ J_ J *A \ 

^(r,t) - p, * p , * ( />, * Pi * J 


* 7T 1 

A A / «•«/ 


J,U (•) , - *• r% -‘^ 

\J* ( f, r«) 


■>/<•( r, r- > I 


The temperature at r = r^ is 


v.,(r„t) = v,(r,,t) 


= r f * . _/ * 

( >/ < fx.pt (. 4> , * Pa)*. 


-XT*/* 


4L • 
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For numerical calculation of (11) and (12), the aerlea In the formulas must 
be computed for sufficiently large n so that the residual; become negligibly small. 
From (11) by putting t = 0 we obtain a relationship 
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p, ( fi. • *») 

* p * 
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J e ( I'N) /(/« f ''i / » ) 







(13) 


Since the lefthand side of (13) is a constant, this relationship can be used to test 
the converge \ce of the series. To evaluate the individual terms of the righthand side 
of (13), must be calculated from (8). Since the coefficient in the equation (8) is 
positive, the roots exist in the intervals where J Q and have opposite signs. 
According to the elementary theory of Bessel function, the signs of J Q and are 
opposite in the interval between and where and are, respectively, the 
n'th positive roots of Jq and Therefore, we have a relationship 


o<xj< J'n^l P n 


(U) 


Since the values of aC and 8 are known and aveilnble in the litereture, the relation- 

n ) n ' 

ship (14) provides an effective way of calculating f by computer. The numerical 
evaluation of (13) showed that, for the ranges of values of the physical constants 
that are encountered in our problem, a summation of first 20 terms leaves the residual 
less than 0.1 % for any value of r that is in the range uir s r^. For the analysis of 
data, y n * s are calculated up to n = 100 and the summation of first 100 terms is used 
to fit the formula to the data. 

The formulas (11) and (12) were obtained for a constant heat flux Fq. We want to 
extend the formulas to the case when the heat flux changes stepwise as a function 0 ** 
time. Let the heat flux be a constant F^ in the time interval between t^ ^ and t^. 

Then, the variation of temperature at time t (> t^) due to the heating is given by 



( 15 ) 




where 


f(r,t) = Vj.fr, t)A 0 
is a function derived from (11). 

The variation of temperature due to the stepwise varying heat flux will be 
obtained by superposing (15) for all time intervals. 

I 


Mr,t) = L v.(r,t) . 

A 4>/ “i — l»“i 

s i £• • . { {< r, A - AO - {<*, t - A ) } 

i* / * • 
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.= r I — — — 

.7/ 1 ♦ *• 


- tx. Z ~^7jr ( Ay:7f7j > (16) 


where t. =t. 

max 
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The temperature at r = is 

? 2 (r 2 ,t) = V^r^t) 
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( 17 ) 


In the case when heat flux F't) is a continuous function of time, the effect of 
heating in the time interval between t^ ^ and t^ can be approximated by that of constant 


* 


i 


heating rate by putting 


rAi 

F(t) dt * F i - (t A - t i-x ) 


Therefore, with an appropriate aelection of time intervals, (16) and (17) will be 
applied to the case when the heat flux varies continuously as a function of time. 
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K. Horn, M. Langveth. Jr. , A Wechsler. J. Winkler, D. Colvin, and S. Keihm k New Technique of 
1 hemal Conductivity Meavuremenl of Lunar Core Sample*. 

A new technique hat been developed to meature thermal conductivity of lunar core tamplct returned 
by Apollo mivtiont. Sm«.e ihcte tamplct arc to be utcd for vanout coentitic ttudict. the mr'turemcnt 
mutt be performed without deteriorating any physical and chemical propcrtic* ot the tamplct. The 
new technique doct not require the extraction of tamplc material from the core tube The lunar 
material in the core tube it mounted in a vertical potmon intide a vacuum cumber and heat dux it 
supplied radially from outvidc of the core lube. The thermal conductivity it determined by analyzing 
the time variation of temperature at the turface of the core lube An experimental vetup hat been 
constructed to tit the Apollo 13 type drive tube core tamplct. Heat it transmuted to the sample by 
radiation from a coaxial cylindrical heater. The heat irantmittion cocthcicnt wit determined ex- 
perimentally for a given temperature difference between the heater and an empty core tube. The 
comparison of the experimentally determined cocthcicnt with the theory of radiative heat transfer 
yielded an estimate of thcrma cmittiviiy c - 0 47 which it reasonable for an anodized aluminum 
turface of the core tube. 1 he solution of the time-dependent heat conduction equation lor a concentric 
composite cylindci (outer cylinder, the core tube, inner cylinder, the sample ot lunar rcgolity) was 
obtained under the conditions that the initial temperature it zero and the influx of heat through the 
outer turface of the vote tube it a continuous function of time. Since the physical properties of the 
core tube are known, a tingle unknown parameter, either thermal conductivity or ditfutivity of lunar 
regolith, can be determined by comparing the theory with the experimental data 

Test experiments have been performed on powdered knippa basalt (grain sizes ranging from 74 to 
149 gz; bulk density. I.39gcm J i and glass beadt (grain size, VOn; bulk density. I.48gcm •*). 
Gaseous pressure in the sample was below 0 78 • 10 " torr for powdered Knippa basalt and below 
1.0 » 10 •* torr for glass beads. For both of these samples, the initial temperature was about 300K, 
which increased about 7k during the measurement The values of thermal conductivity obtained by 
the present method were 1. 10 ■ 10 - cal cm -1 s 1 k for po edered knippa basalt and 0 41 10 ‘cal 

cm* 1 s *k for glass beads These values are consistent with the values obtained by other standard 
techniques cn these mxtenals. It is considered that the new technique is particularly suitable for the 
measurement of thermal conductivity of lunar core samples for the following reasons: III The 
measurement is mace without disturbing the sample mechanically. It is not necessary to dissect the 
core tube and extract the lunar maicnal for the experiment (2) The sample is heated in an indirect 
manner. Heat is transmitted to the core tube by radiation to avoid a direct contact to the heat source 
The sample temperature need not be raised more than 10k to obtain reliable data (3i 1 he attachment 
of a temperature sensor to the outer suilace of the core tube is the only mechanical contact required 
for the meavuremenl. This will greatly reduce possible chances of chemical contamination to the 
sample during the measuremetr. 

This work was supported by NASA gram NGR 33-008-169. 



